Quantum radiation force on a moving mirror with Dirichlet and Neumann boundary 
conditions at vacuum, finite temperature and coherent states. 
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We consider a real massless scalar field in a two-dimensional spacetime, satisfying Dirichlet or 
Neumann boundary condition at the instantaneous position of a moving boundary. For a relativistic 
law of motion, we show that Dirichlet and Neumann boundary conditions yield the same radiation 
force on a moving mirror when the initial field state is invariant under time translations. We obtain 
the exact formulas for the energy density of the field and the radiation force on the boundary for 
~qq ' vacuum, thermal and coherent state. In the nonrelativistic limit, our results coincide with those 

' found in the literature. 

o '■ 

In the 70s decade the first works investigating the quantum problem of the radiation emitted by moving mirrors in 
vacuum werepublished, motivated by the investigation of particle creation in the nonstationary universe (see Refs. 
iS ! EE 0, 1 1 H S ) ■ Fullling and Davies 1 studied the moving mirror radiation problem in the context of a real scalar 
■ field in a two dimensional spacetime. They obtained an exact formula for the finite physical part of the expected value 
of the energy-momentum tensor, assuming the initial state as the vacuum. Their results revealed that the radiation is 
originated at the mirror and propagates away from it. Ford and Vilenkin Q developed a perturbative method which 
can be applied to mirrors moving in small displacements and with nonrelativistic velocities. In this approximation, 
they obtained, for a real scalar field in a two dimensional spacetime, that the radiation force is proportional to the 
third time derivative of the mirror law of motion, which is the nonrelativistic limit of the result obtained in Ref. 
[H]. Ford and Vilenkin also applied their method to a scalar field in four-dimensional spacetime, obtaining a force 
proportional to the fifth time derivative of the displacement of the mirror. Davies and Fulling [j| (see also Ref. @) 
found a class of hyperbolic trajectories of a mirror moving in a 1+1 Minkowski spacetime, for which the emitted 
radiation corresponds to a thermal spectrum, that could be registered by an inertial detector. Moreover, several 
works have focused on the problem of moving mirrors placed in a thermal bath (considered as the "in" state, instead 
of the vacuum) . Jaekel and Reynaud obtained for the scalar field in 1 + 1 dimensions the thermal contribution 
to the dissipative force proportional to the velocity of the mirror, valid in the nonrelativistic limit. Thermal effects 
C*~) , have also been considered in Refs. [Ill, E2, E!> EH EE EH • The coherent state is another initial field state which has 
been considere d HtL Il8j , as well as the superposition of coherent sates used to study decoherence via the dynamical 
Casimir effect 19]. Furthermore, attention has been given to the role that boundary conditions play on the dynamical 
£f~) ■ Casimir effect. In the static Casimir effect, different boundary conditions can change the sign of the Casimir force 



20]. Applications of this change, for instance, in the building of nanoelectromechanical systems have been discussed 
21 j . The role of Dirichlet and Neumann conditions on the static Casimir force has been investigated in Refs. [22[ in 
the context of the scalar field. Recently, several works have investigated the influence of different boundary conditions 
on the dynamical Casimir effect [13, [23[. In this context, Farina, Maia Neto and one of the present authors, showed 
that Dirichlet and Neumann boundary conditions yield the same force on a moving mirror when the initial field state 
is symmetrical under time translations jl7j |. However, the validity of this conclusion, obtained in the context of the 
Ford- Vilenkin approach Q, was restricted to the nonrelativistic and small displacement approximations. 

The exact method of Fulling-Davies Q is frequently considered as less flexible than the perturbative methods 
because it is based on conformal transformations and applicable only for the two-dimensional spacetime [3, EH, 
whereas the perturbative approach can be extended to higher dimensions. On the other hand, the Fulling-Davies 
method enables to obtain exact results for the scalar field in the two-dimensional spacetime. In the present work, 
instead of following the approximate approach considered in Ref. [TtJ , we use the exact approach proposed by Fulling 
and Davies [H, and show that, also for relativistic laws of motion, Dirichlet and Neumann boundary conditions yield 
the same radiation force on a moving mirror when the initial field state is symmetrical under time translations. In this 
context, we calculate for a real massless scalar field in a two-dimensional spacetime the exact (relativistic) formulas 
for the energy density of the field and for the force on the mirror, for Dirichlet and Neumann conditions, at vacuum, 
finite temperature and coherent initial states of the field. In the nonrelativistic limit, these exact formulas lead to the 
approximate results found in Refs. [ToL fl7j ■ 

We start by reviewing the exact field solution for Dirichlet and Neumann dynamical boundary conditions 0, 0, 
0]. Let us consider the field satisfying the Klein-Gordon equation (we assume throughout this paper H = c = 
1): (<9 t 2 - dl) cj) (t, x) = 0, and obeying Dirichlet (<£'(*', x')\ boundarv = 0) or Neumann (d x ><f)' (f ,x')\ bound ary = ) 
condition, taken in the instantaneously co-moving Lorentz frame, at the moving boundary position x — z(t). We 
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FIG. 1: Moving mirror trajectory. The dashed lines are null-lines separating the regions I from II, and III from IV. 



examine the particular set of mirror trajectories for which z(t < 0) = 0, as shown in Fig. [T] Using the appropriate 
Lorentz transformation, Dirichlet and Neumann conditions can be written in terms of quantities in the laboratory 
inertial frame as follows: <j>[t, z(t)] = and {[i(t)d t + d x ](f>(t,x)}\ x=z(t - ) = 0, respectively. The field solution can be 
obtained by exploiting the conformal invariance of the Klein-Gordon equation. Considering the conformal coordinate 
transformations t — x = f (w — s) and t + x = g (w + s) , the scalar wave equation in 1+1 dimensions remains invariant: 
(<9^ — df) <p (w, s) = 0. Functions / and g can be set up, for a large class of laws of motion z(t), so that the curve 
s — coincides with the boundary trajectory x — z(t): [t,z(t)] — > {w,0) Then: <f>[t,z(t)] = — > 4>(w,0) = 
and {[z(t)dt + d x ] 4> (t, x)}\ x=z ^ = — » [d s 4>(w, s)] s =o = 0. The mode solution of the wave equation with static Dirichlet 
or Neumann boundary conditions in (to, s) space are well known so that, coming back to (t, x) coordinates, we get: 
(j) (i, x) — Jq 00 did [a u (j)u + aj,0*] , where 



4>u (t,x) = (47TW)" 



7e 



-iw!rr(v) _j_ * — iuJ7rp(u) 



(1) 



form a complete set of positive-frequency solutions, and u = t — x, v = t + x. In Eq. ([T]), we introduce a notation 
which enables us to put into a single formula the solutions for Dirichlet and Neumann boundary conditions, and also 
the solutions for the right hand side (regions I and II in Fig[T]) and the left side (regions III and IV) of the mirror. 
In this sense, for 7=1, Eq. (fTJ) gives the Neumann solution, whereas for 7 = i we have the solution for Dirichlet 
boundary condition. For the regions I and II showed in Fig. [TJ r (v) = v and 2t(u) —11 = f^ 1 (u) = p (u), where t(u) 
can be obtained from t(u) — z[r(u)] = u; for the regions III and IV: p (u) = u and 2r(u) — v = g^ 1 (u) = r (y), where 
t(v) + z[t(v)] = v. As causality requires, the held in the regions I and IV, represented by cj>o, is not affected by the 
boundary motion Q, so that p and r are also chosen to be identity functions in these static regions. Hereafter we 
consider the averages (...) taken over an arbitrary initial field state (regions I and IV) assumed here, for simplicity, 
as being the same one for both sides of the mirror. We continue our analysis writing the correlator function C = 
{4>o(t, x)(f>o(t', x')) split in the following manner: C = C vac +C/ata) +C(aa), where 



C 



du}Fi(u,u', |7|,7)| a 



OO pOQ 



C{aia)=J J dojdu' (al ) ,ajjj 7 -i(ui,u)\\j\ ) j) + c.c. 



C(aa) — 



00 /-OO 



^0 



and 



dujduj' {d^'du) Fi(uj, —u',7, 17I) + ' 



^(w.w'.ftA) = e- l ^ t -"' t ')/(4nV^?)\p 2 e- l ^ x - UJ ' x '') 
+\ 2 e - % ( u>xWx ') + c.c. . 

From these equations, we see that, in the presence of the boundaries, C va c,C^a) an d are not symmetric under 
space translations. The term C vac is symmetric under time translation, whereas C^ta> is symmetric if /a^,a w \ cx 6(u>'— 
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uS). On the other hand C(aa), if non-null, is not symmetric under time translations. In the context of the perturbative 
approach Q , the radiation force can be given in terms of correlation functions depending on the unperturbed field 
operador <j> : \d x d x iC] T ,_ T _n and [C] x ,_ x _ for Dirichlet and Neumann boundary conditions respectively [1,0]. As 
shown in Ref. [17] . the parts of the force related to C vac and C^ta) are the same for Dirichlet and Neumann boundary 
conditions. On the other hand, the part related to C/&&\ is different by a sign. In other words, in Ref. [17j it is shown 
that the difference between the force acting on Dirichlet and Neumann emerges, in the context of the non-relativistic 
mirror motion with small amplitude, from the part of the non-perturbed correlator C which is non-invariant under 
time translations. 

Now we will investigate this problem in the context of an exact approach, starting from the field solution fl} and 
calculating the exact formulas for the expected value of the energy density operator T — (Tbo(t, x)), and the net force 
F (t) acting on the moving boundary defined by (since Too = ^li m this model) : F (t) = T [t, z(t)] — 1 [t, z(t)f + \ 
where the superscript "+" indicates the regions I and II, whereas "-" indicates the regions III and IV in Fig. [T] Let 
us split T, writing T = % ac + T^a) + ?(&&)■, where: 



/>oo 

% ac = 1/2 / duj T%{w,u)\ 
Jo 



71 



and 



/>oo />oo 

Jo Jo 

/>00 r-OG 

T~{aa)=-l/2 / du;duj'(a U j>a ul )J 7 2(u 7 uj',j,"{*) + c.c., 
Jo Jo 

,p,X) = — — \p [r (y)\ e ( >( > 
+ X 2 [p 1 {u)f e- l{w ' w ' )p{u) Y 

The first term T vac is the local energy density related to the vacuum state, which is divergent. After regularization 
and renormalization (see Ref. Q), T vac can be redefined as the renormalized local energy density: 

% ac = H7I7247T [p'" (u)/p' (u) - (3/2)p" {uf /p> (uf 
+p(u) -> r(v), 

where the object appearing in the brackets is known as the Schwarzian derivative. We remark that, in our notation, 
for the right side of the mirror we have r(v) — v, so that this equation gives back the formula found in Ref. @ for 
Dirichlet, and in addition we verify that the same formula is valid for the Neumann boundary condition. Now we 
investigate the net radiation force F acting on the moving mirror. Let us consider F = F vac + iTata) + F(aa)i where: 

Fvac = | 7 | 2 (l + i 2 ) 

x (z 2 i/2^)/(l - i 2 ) 4 + («/67r)/(l - i 2 ) 3 

poo poo 

F/xia\ = 2 



F{aa)=-l / dudw' (a^a^) T z (u),u' : 7*, 7) + C.C. , 



Jo 



and 



47T ^ (1 - z) 

+p 2 e -*(w-«')[t-«(*)] 

- [(z,z,p,p) -> (-z, -z,r,\)]} . 



4 




FIG. 2: Normalized exact thermal force \F^\/ot (solid line) and approximate thermal force \F^?\/(Tt (dashed line), both 
valid for Dirichlet and Neumann boundary conditions, plotted as functions of the mirror velocity z. 

We see that F vac and Fi^gs depend on I7I 2 , which has the same value for Dirichlet and Neumann conditions. On 
the other hand F/ss) depends on 7* 2 ( or 7 2 ), which differs by a sign in Dirichlet and Neumann cases. Noting 
that the term C^a) 7^ 0) is non-symmetric under time translations, we can generalize the result found in Ref. 

[]~7l ] . concluding that for a general (relativistic) law of motion, Dirichlet and Neumann boundary conditions yield the 
same radiation force on a moving mirror when the initial field state is symmetric under time translations. In the 
non-relativistic limit, we recover F vac (t) « |7| 2 z'/67r, found in Ref. [8j for Dirichlet, and in Ref. 17] for Neumann 
condition. 

Let us examine the radiation force on moving mirrors when there are real particles in the initial state of the 
field. We start with the thermal bath with temperature T, which is an example of invariant field state under time 

translations. For this state we need to take into account that (d$.,aj\ = ¥i{lo)5 (u> — uo 1 ) where n{u>) = l/{e hu / T — 1), 

(T) 
(a+a)' 



with the Boltzman constant equal to 1. We get T^a) = and T^a): renamed as the energy density T^iy is given 
by: 7#> = 



7| 2 7rT 2 /12 r' (v) +p' (u) . The force F^a) — 0, whereas -PVta)) renamed as the net thermal force 



F^ T \ is given by: 



F (T) = _ 



, (1 + ^ 2 ) 
(1-i 2 ) 2 



Dz 



2n+l 



where ar = 2|7| 2 7rT 2 /3 is the viscosity coefficient. This exact formula is a generalization of the one obtained in Ref. 
(To| . For non-relativistic velocities related, for instance, to mechanical motions of the mirror (i ~ 10~ 8 ) or to the 
simulation of the mirror motion by changing the reflectivity of a semiconductor by irradiation from laser (i ~ 10 -3 ) 

[24]], the series can be truncated in n — 0, leading to the approximate formula: F^ w = —2\j\ 2 ttT 2 z/3, 

in agreement with Ref. [l(| (for Dirichlet), and also with Ref. [13] (for Neumann boundary condition). In the 
nonrelativistic context, this approximate formula is in good agreement with the exact value. For relativistically moving 
mirrors (for instance, relativistic electron beam may be an embodiment of a relativistic mirror [HI), corrections to 
the approximate formula can become necessary, as it can be seen in Fig. [21 

Let us now consider the coherent state, as an example of a non invariant state under time translations. The 
coherent state of amplitude a is defined as an eigenstate of the annihilation operator: \a) = a5 (w — loq) \a) , where 
a = \a\ exp(i#) and ujq is the frequency of the excited mode. For this case, we have 7^t&) an d T^a) relabeled as Xgt a \ 



and respectively, and given by: 



% 



(a) 
dtd) 



7 | 2 /(27rH{H 2 [r' (vf + p' (u) 2 ]} 



(aa) ^ 1 



{a 2 [ 7 2 \r' {v)f e-^W 
+7* 2 [p' (u)} 2 e- 2luJop{u) 



+ c.c 



}• 
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FIG. 3: We plot the exact coherent force as function of time, for the Dirichlet case. We consider a = exp (in/2), ujq — 10, and 
the following values of the velocity: — 1CP 8 (dotted line), — 10~ 2 (dashed line) and —8 x 1CP 1 (solid line). The dotted and solid 
lines exhibit the coherent force multiplied by the factors 5 x 10 and 1/125 respectively. 



The exact forces i^ata) and F^a), relabeled as the coherent forces -F^ta) ana - -^(aa) respectively, are given by: 



F, 



ata) = —"0 \a\ z(l + z)/(l-z) , 



7 («) 

(dd) 



.^o . |a e - ai(wo t- 9 ) 

47T 



x S 7 



e 2iuj z(t) I 1 _ e -2iu a z(t) 



1 + Z 

1-z 



1 + Z 



e —2itj z(t) _ e 2iw z(t) 



C.C. 



If we consider simultaneously nonrelativistic velocities and small displacements (in the sense considered in Rcf. 17]), 
according to what is required by Ford-Vilenkin approach Q, the force F^ — F^ &) + ^/aa) can ^ e approximated as: 



FW » -^ H 2 {2| 7 | 2 i(0-(7 2 + 7* 2 ) 

x [cos (2w i - 20) z (t) - sin (2w t - 26») w z (*)]} ■ 

In Fig. [3] we plot the exact coherent force as a function of time for the Dirichlet boundary condition and different 
values of the mirror velocity. Assuming the mirror moving with uniform velocity toward the negative direction of the 
x-axis, the Fig. [3] shows the force oscillating F^ and the graph shifting to the positive region of the vertical axis as 
the mirror velocity grows, becoming the force more intense and opposite to the motion. For the Neumann boundary 
condition, the force oscillates in a different manner, but exhibits analogous shift for relativistic velocities. 

In summary, focusing on the advantages of the Fulling-Davies approach for the case of a massless scalar field in 1 + 1 
dimensions, in the present paper we showed the exact dynamical Casimir force acting on a moving boundary under 
Neumann condition, with the vacuum as the initial state, generalizing the non- relativistic result found in Ref. [13] • 
For the thermal field, considering both Dirichlet and Neumann conditions, we wrote the exact formula for the thermal 
force, generalizing the approximate formula found in Ref. [l(| and also, for instance, in Refs. [HI. [l2l Il5l [1(1 fl7j . For 
the coherent initial state, we found exact formulas for the radiation force, which are different if we consider Dirichlet 
or Neumann condition, generalizing the perturbative formulas found in Ref. [13] • Finally, we extended to a general 
(relativistic) law of motion the conclusion found in the literature [17| that Dirichlet and Neumann boundary conditions 
yield the same radiation force on a moving mirror when the initial field state is invariant under time translations. 

We acknowledge L. C. B. Crispino and P. A. Maia Neto for careful reading of this paper. This work was supported 
by CNPq - Brazil. 
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